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Abstract 

We employ the analytic Bethe Anzats to construct eigenvalues of transfer matrices with 
finite-dimensional atypical representations in the auxiliary space for the putative long-range spin 
chain encoding anomalous dimensions of all composite single-trace gauge invariant operators of 
the maximally supersymmetric Yang-Mills theory. They obey an infinite fusion hierarchy which 
can be reduced to a finite set of integral relations for a minimal set of transfer matrices. This 
set is used to derive a finite systems of functional equations for eigenvalues of nested Baxter 
polynomials. 



1 Introduction 



Up until very recent years, only one- dimensional quantum field theories entertained nonpertur- 
bative analytical considerations in coupling constant [TJ [2]. The reason behind availability of 
these results is largely due to exact solubility of corresponding models [HI HI [2] . Four-dimensional 
non-supersymmetric gauge theories were also known to possess integrable structures. In certain 
limits, their dynamics effectively becomes low- dimensional and one observes appearance of struc- 
tures similar to the interacting spin chains. The high-energy limit of QCD scattering is described 
by a two-dimensional reggeon theory [5] , whose quantum mechanical Hamiltonian was shown to 
possess holomorphic separability [6] with each component being identified with the Hamiltonian 
of the one-dimensional noncompact sl(2, C) Heisenberg magnet 0, [7]. In a different limit, the 
one-loop renormalization group evolution of quasipartonic Wilson operators [8] is governed by a 
pairwise light-cone dilatation operator. In the multi-color limit, only nearest-neighbor interac- 
tions survive in it and its maximal-helicity subsector is mapped to the XXX s/(2,M) spin chain 
O [TOj HH [12] . In both cases, the models were solved by powerful techniques such as Bethe Ansatz 
[T31 14] and Baxter Q— operator [T4] . 

The integrable spin chain underlying the QCD maximal-helicity dilatation operator gets ex- 
tended to sectors with other field contents in gauge theories enjoying supersymmetry [T5J [161 
[T7J HH HU [201 ■ In the maximally supersymmetric Yang-Mills theory, anomalous dimensions 
of all Wilson operators [T5J, [THl [2T] are encoded in the underlying su(2, 2|4) spin chain [19] that 
can be diagonalized by means of the nested Bethe Ansatz [22]. The symmetry algebra of the 
magnet reflects the classical symmetry of the gauge theory. Moreover, the Af = 4 supersym- 
metric quantum field theory, being anomaly-free, remains superconformal even at the quantum 
level such that the absence of bound states in its spectrum warrants anomalous dimensions to 
be the relevant physical observables. Thus their calculation implies solution of the field theory 
in question. What is even more extraordinary is that in the M = 4 superYang-Mills theory, 
the one-loop integrability, manifesting itself through spectral degeneracies, persists beyond lead- 
ing orders of perturbation theory. This was demonstrated by explicit multi-loop diagrammatic 
calculations in compact [17] and noncompact [231 121] subsectors of the theory. This evidence, 
together with multiloop field-theoretical calculations for twist-two operators [251 [26] , pointed out 
that anomalous dimensions of corresponding gauge invariant operators can be deduced from a 
kind of a deformed Bethe Ansatz as was shown in Ref. [27]. Together with insights [28, 29] into 
exact solvability of this gauge theory coming from its conjectured duality to the type IIB string 
theory on the AdSsxS 5 background [30], a set of all-order nested Bethe Ansatz equations was 
suggested in Ref. [31]. 

The conjectured long-range Bethe Ansatz equations are exact only for infinitely long spin 
chains. For finite-size systems they become invalid when the interaction range becomes longer 
than the number of fields in the gauge-invariant operator. Thus they are only asymptotic in 
character (see, e.g., Ref. [32]) and this calls for consistent incorporation of finite-size effects into 
the current treatment. One approach to resolve this problem is known as the Thermodynamics 
Bethe Ansatz (TBA) [331 EI] extended to lattice systems in Refs. [35, 36J. A pivotal role in 
the derivation of TBA equations, which is independent of the notorious string hypothesis, is 
played by the functional algebraic relations among the hierarchy of fused row-to-row higher spin 
transfer matrices [371 ESI EH HD] • For simplest cases, these relations were shown to have a deep 
mathematical origin in exact sequence of Yangian modules [39] . The resulting fusion hierarchy 
yields an infinite system of coupled equations which can be consistently truncated only for finite 
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Figure 1: Kac-Dynkin diagram used in the construction of transfer matrices. 



temperature and is hard to practically use otherwise for vanishing temperatures |41j . However, it 
was demonstrated by Takahashi [32], and extended to superspin chains by Tsuboi [JT], that the 
infinite hierarchy of TBA equations can be reduced to a finite set of integral equations for a finite 
number of fused transfer matrices. This represent an alternative and more practical approach to 
study the thermodynamics of integrable models. 

In this paper we address the question of constructing fusion hierarchies among transfer matri- 
ces for the putative su(2, 2|4) long-range spin chain of the maximally supersymmetric Yang-Mills 
theory and their use for the solution of the spectral problem. The peculiarities of the model 
exposed in renormalization of the spectral parameter and various dressing factors make the an- 
alytical properties of fused transfer matrices quite complicated to allow for generalization of 
Takahashi equations. Consequently, we focus on the derivation of closed systems of Baxter 
equations involving minimal sets of fused transfer matrices. These can be diagonalized either 
order-by-order of perturbation theory [131 EH SSI HE] or used for strong coupling analyses as was 
demonstrated earlier for the long-range rank-one sl(2) spin chain [4"T] . 



2 Analytic Bethe Ansatz 

The main ingredient of the current framework are row-to-row transfer matrices that encode full 
set of mutually commuting conserved charges of the chain. For short-range spin chains, the 
R— matrix formalism provides the relevant framework for systematic construction of commuting 
transfer matrices. We possess however only the conjectured nested Bethe Ansatz equations [31] 
for the underlying long-range magnet. The lack of a systematic microscopic approach does not 
prevent one from finding the eigenvalues of fused transfer matrices though. The main role in 
this situation is played by the analytic Bethe Ansatz [H] which relies on general macroscopic 
properties of the system such as analyticity. 

For the spin chain based on the high-rank su(2, 2|4) superalgebra, there are several equivalent 
ways to formulate the nested Bethe Ansatz [19J [50] . They reflect different choices of its simple 
root system {ck p |p = 1, . . ., 7} [51]. For the rational short-range R— matrix, which generates the 
spectrum of one-loop anomalous dimensions, they take the following generic form 

',>) _ in n \ L 7 n " ii (p) - iS q) 4- iA 

u 0,k 2 W P \ TT TT u 0,k u 0,j ^ T^PQ (r) i \ 

+ ~ 1111 „0») _„(«)_ i A ' 1 ' 

u 0,k'2 W P/ q=lj=l u 0,k u 0,j 2 PQ 

with the structure of the algebra encoded in the Cartan matrix A pq = (ct p \ct q ) and the representa- 
tion on the spin-chain sites determined by the Kac-Dynkin labels, w p = (a p \fi) expressed by the 
corresponding weight vector fi. However, out of many equivalent Kac-Dynkin diagrams only four 
allow for an all-order generalization [31], i.e., consistent deformation of ( 12.11) with the 't Hooft 
coupling g = g YM ^/W c /{2T:). They correspond to the ground state \fl) built from the complex 
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scalar Z— fields of the theory, \Q) = tr[Z L (0)], which is protected from quantum corrections to 
all order of perturbation theory. We choose to develop our formalism based on the Kac-Dynkin 
diagram displayed in Fig. [T] with its central node corresponding to the light-cone si (2) subal- 
gebra formed by projections of the generators of the Lorentz transformations, translations and 
conformal boosts, 

iWT + = L 1 1 + L i i , ity + = V2F li , iST = \/2K il , (2.2) 

respectively. The pseudovacuum of the second level Bethe Ansatz corresponds the derivative 
sector |0') = ti[D^Z L ] universal to all gauge theories. The only nontrivial Kac-Dynkin label 
of the gauge-theory spin chain is associated with the central node, W4 = —1. The complete 
Serre-Chevalley basis is spelled out in Appendix lAl 

Transfer matrices are supertraces of monodromy matrices with a chosen representation of 
the symmetry algebra in the auxiliary space. Thus their eigenvalue formulas are expected to 
be a sum of terms corresponding to each components of a given representation labelled by the 
Young diagram as was demonstrated in Refs. [37J EH] and [40J for classical and graded alge- 
bras, respectively. Thus the construction is based on the identification of the elementary Young 
supertableaux which depend on the spectral parameter and read for the long-range su(2, 2|4) 
magnet 




Here each box corresponds to the respective bond connecting the nodes of the graph in Fig. [T] 
with assigned bosonic 1=4 = 5 = 8 = and fermionic 2 = 3 = 6 = 7=1 gradings. These 

1 See Appendix [El for Young tableaux corresponding to the Kac-Dynkin diagram with the FBBFFBBF grading, 
i.e., with the central compact sit (2) node. 
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tableaux are built from bare and dressed Baxter polynomials, 

n p n p 

gw(u) = n (« - 4 p) ) , ^ (p) («) = n (* - . (2-5) 

k=l k=l 

expressed in terms of, respectively, bare u and renormalized spectral parameter x = x[u] = 
\{u + \Ju 2 — g 2 ) [29], and assumed shorthand notations = u ± | and x 1 * 1 = x^]. The 

Bethe roots admit an infinite series expansion in the 't Hooft coupling, = uf\g) = 

u ok + 9 2vj{ fk + ®(9 4 )- The exponential factors correspond to dressing factors which emerge at 
higher orders of perturbation theory and read 



f l dt lnQ^^-gt) ( 

Li \ u + gt J ' 1 ' } 



where rj takes three possible values ±, 0. For rj — 0, it can be represented as a ratio of the bare 
and dressed Baxter polynomials 

VV) =ln g!M. (2.7) 
2 v ' Q(p)(«±) V 7 

For the momentum carrying Bethe roots, there is an additional contribution Q(u) 

A±(w) =<ri 4) (w) -©(«), (2.8) 

known as the magnon dressing phase which is responsible for smooth interpolation between the 
weak- and strong-coupling expansions [52"| 1531 15^] and reads in terms of the Baxter polynomials 



x / — * lnfl+ ^ / T V (2.9) 

■/qMoo] 27U sinlr (vr/t) V 4xx i K + # S J / V Axx l K ~ 9s\ J 

Introducing a symbolic notation for the differential operator involving the single-box Young 
supertableau introduced above, labelled by the flavor index a and the spectral parameter u, 



V aA±] = l±\a\ u e- l9 \ (2.10) 



one can write the generating functions [37J EH EHJ HO] for the eigenvalues of transfer matrices with 
atypical finite-dimensional antisymmetric (l a ) and symmetric (s) representations in the auxiliary 

space, 

oo 

^ l[+ ]^f +] ^ 6 -f +] P 5 , [+ ]^4j + ]^M^ [ 1 + ]^i,[+] = E *M (^ [1_al ) e " m9 " > (2- 11 ) 

a=0 

oo 



s=0 



where 

u [±n] = u ± in , x [±n] = x[w [±n] ] , (2.13) 

which are generalizations of the previously introduced conventions w^ 1 ' = u ± . In Eqs. (12.111) and 
(I2.12p . the powers p a = 1 — 2a of the multiplicative factors in the left-hand side of the equations 
reflect the grading of the bonds of the Kac-Dynkin diagram. The residues of apparent poles 
positioned at zeroes of the Baxter polynomials in these transfer matrices vanish upon the use of 
the long-range nested Bethe Ansatz equations of Ref . [31] . 

The momentum carrying Baxter polynomial Q^ 4 \u) determines the anomalous dimensions 

of local Wilson operators with quantum numbers defined by the labels of the Kac-Dynkin dia- 
gram, Fig. [H as demonstrated in Appendix |A] Solutions corresponding to the anomalous dimen- 
sions of Wilson operators are selected by the vanishing quasimomentum 6 condition imposed on 
the states of the chain 

„ if 1 dt , Q (4) (+f-^) , s 

i0 = - I In 2 — ^4 = 0. 2.15 

For the sm(2, 2|4) magnet, the strategy to determine the eigenvalues Q^'(u) consists in se- 
lecting the first seven fused transfer matrices, either t\ a ](u) or P- S ^(u), and subsequent solution 
of the resulting system of equations for the zeroes of transfer matrices and Baxter polynomials. 
For the Kac-Dynkin diagram we have chosen for our analysis, one immediately finds that t\ a ](u) 
possess kinematical zeroes 

a- 2 

t [a] (u) = H(J a - 2 ~^) L T [a] (u), (2.16) 

3=0 

contrary to the ones with symmetric auxiliary space t^(u). Therefore, while for the former one 
has sufficient number of equations to determine all unknowns, it is no longer the case for the 
symmetric transfer matrices that contain a redundant set of conserved charges. The latter are 
not completely independent though and can be related by means of nonlinear integral equations 
obeyed by t^(u) [121 HI], as we demonstrate at one loop in Section HI We will study the 
diagonalization procedure below on examples of low-rank subsectors. It will turn out that due 
to simplicity of corresponding fused transfer matrices, the rank-two s/(2|l) subsector allows for 
complete resolution of nesting. Above it, the nesting is intrinsic, however, it may disappear in 
specific cases once the consideration is restricted to the lowest order of perturbation theory. 



3 T-system 

The antisymmetric and symmetric transfer matrices (12. lip and (I2.12p discussed in the previous 
Section obey an infinite set of algebraic functional relations. To find them, we introduce at 
first transfer matrices with the auxiliary space labelled by a general Young superdiagranj^ Y(s) 
defined by the partitioning s = {s 1; s 2 , . . . } with s\ > s 2 > ■ ■ ■ [37]. To this end, one labels 

2 Our analysis can be extended to skew Young supertableaux as was done in Ref. [46] for the sZ (2 1 2) sector, 
however, since these will not be relevant for our subsequent considerations and they will not be presently addressed. 
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all boxes of the Young superdiagram Y(s), starting with the one in the top left corner, with 
a pair of integers (j,k), corresponding to rows and columns, respectively. Then one defines a 
set of admissible Young supertableaux Y a (s) by assigning one of eight a(j,k) "flavor" indices 
to each box of the diagram Y(s) and distributing them according to the ordering conditions: 
a(j, k) < a(j,k + l) and a(j, k) < a(j + 1, k). These strong ordering conditions turn into weaker 
inequalities when the indices on adjacent boxes have coincident gradings, namely, for 

• bosonic grading a = 0: 

a(j,k) <a(j,k + l), a(j,k) < a(j + l,k) ; (3.1) 

• fermionic grading a = 1: 

a(j, k) < a(j, k + 1) , a(j, k) < a{j + 1, k) . (3.2) 

The transfer matrix for a rectangular s = {s a } Young superdiagram Y(s) can be easily con- 
structed according to this set of rules and reads 



Y a a(j,k)£Y a 



a(j,k) 



,[s— a+2j — 2k] 



(3.3) 



Obviously, by construction, tP?(u) = t[ a ](u) and tj^(u) = t^(u). This general rectangular 

transfer matrix t$(u) admits the Bazhanov-Reshetikhin determinant representation [33 SO] 
(see also Ref. [56] for a recent discussion) in their terms, 

t\$(u) = Jet t [a ^ +k] (u^ +1 -^) = Jet if+i-*}(„li+*— i]) , (3.4) 

with boundary conditions v- n *(u) = t\ n ](u) = for n < 0. Making use of the Desnanot-Jacobi 
identity for determinants, one finds the T— system of equations [351 ESI [551 HO] (see also [57] ) 
obeyed by rectangular transfer matrices, 

4V)i s ]V) = t\^\u)ti:j\u) + t&wtf^cu) , (3.5) 

with the following truncation properties which will be crucial for the analysis in the following 
Section 

tj>>( u ) = o ; f or a , s > TV + 1 , (3.6) 

and the duality relational 

t { m n} (u)=t { ^ +n] (u), (3.7) 
with TV = 4 for su(2, 2|4) and n being non-negative integer. 



3 We kept AT arbitrary so that these results could be applicable to other sl(Af\J\f) algebras. 
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4 Takahashi integral equations 



As we pointed out at the end of Section [2, the system of seven equations for the Baxter polynomi- 
als Q( p \u) involving the symmetric transfer matrices t* a *(u) yields an underdetermined system of 
equations for their zeroes. To reduce the number of unknowns, they have to be supplemented by 
relations between redundant conserved charges. The starting point in the derivation of nonlinear 
equations, the so-called Takahashi equations, for the symmetric transfer matrices is the fusion 
relation 

H (x^- 1 - 2 ^ (t^(u + )t^(u-) -t^iu^Hu)) = T^ } (u) , (4.1) 

j=0 



where we have used the following factorization property of the fused transfer matrix tj^(u) in 
its right-hand side 

U(x lMj r L t$(u)=T$(u), (4.2) 

3=0 

which enables the very possibility to derive nonlinear equations. An obvious way to relate the 
charges entering the expansion of transfer matrices is by using the fact that the renormalized 
transfer matrix T^){u) is a polynomial in the renormalized spectral parameter, so that the 
residues of poles in the left-hand side have to vanish on their own order by order in perturbation 
theory. At leading order, this can be summarized in a closed set of nonlinear integral equations 
as we demonstrate below. 

One notices that the fusion hierarchy (14.11) can be rewritten as 

The crucial observation is that the left-hand side of this equation possesses poles (at most of order 
L) at u — — ~(s — 1 — 2j) with j = 0, . . . , s — 1, while the right-hand side does not. Expanding 
the first term in the left-hand side as a sum of poles, 

n (^-^r L 4 s V) = ei>(*) ( ulMj] v k > (4.4) 

j=0 j=0 k=l 

one finds the coefficients aj(k) by integrating both sides of this equation in the vicinity of sin- 
gularities. Here we used the fact that the coefficient of the leading asymptotic term at large u 
vanishes, t^(u) = ■ u sL + (9(u sL ~ 1 ). Performing the summation over k, one deduces 

ne^-MTiJ'V) = / ^^il^-'-'T^V), (4.5) 

where the integration contour C s is shown in Fig. [2j Substituting ( 14.31) into the above equation, 
we can neglect its right-hand side since it does not have poles at u — — |(s — 1 — 2j) and obtain 
a system of integral equation for eigenvalues of transfer matrices 

1 ; Jc.2*iu-vfl \vl-i-*l) t\ s} (v-) 1 ^ 
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Figure 2: The integration contour in Eq. (14.51) . 



As it stands, this is an infinite set of nonlinear integral equations. However, as a consequence of 
the identity (13.71) . it truncates. Namely, using the Bazhanov-Reshetikhin formula for both sides 
of this equality, we can solve it for tg 2A ^ {u) and find the latter in terms of the transfer matrices 
t { s} {u) with 1 < s < 2N, 



det (1 - 4'-^\uU^-^) 

l<7,fc<JV 



det (1 

l<7,fc<AAfl 



(4.7) 



det $- k+ * r} (uV+ k -W)+ det 4 j - k+X - 1} (^+ k -^) 

l<j,k<Af-l U l<j,k<Af U 



As} 



Such that one gets a finite system of integral equations involving only transfer matrices 
1 < s < 2Af — 1 with M = 4. Having reduced the number of conserved charges, the system of 
equations for eigenvalues Q^(u) can be solved. 



5 Baxter equations 

In this Section we derive closed systems of TQ-relations, known as Baxter equations, for Baxter 
polynomials in low-rank subsectors of the su(2, 2|4) long-range spin chain. Contrary to previous 
studies [121 HH] instead of dealing with conjugate transfer matrices, we will use fused transfer 
matrices with increasing spin in the auxiliary space. 

First of all, removing all excitations except the ones associated with the central noncompact 
node of the Dynkin diagram, we immediately find that the transfer matrix t^(u) = t^(u) yields 
the known Baxter equation of the long-range sl(2) magnet (131 US]. The ground state of this 
chain is built from the product of complex scalar fields corresponding to Z = 43 in conventions 
of Ref. [58], and the light-cone covariant derivatives as elementary excitations propagating on it. 
The iterative numerical solution order-by-order in perturbation theory is given in Appendix [C] 
Let us proceed to higher rank sectors. 
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5.1 s/(2|l) sector 

Extending considerations to the case of nonvanishing excitations on the first fermionic node 
adjacent to the central one, either left or right, we obtain the closed s/(2|l) sector. Now in 
addition to covariant derivatives, it possesses fermionic excitations, which correspond to either 
helicity up A 1 or helicity down A4 gauginos, depending on the choice of the fermionic node. Let 
us address the case with the nontrivial nodes three and four, thus Q^w) = 1 for p 7^ 3,4. 

One can derive Baxter equations either using symmetric or antisymmetric transfer matri- 
ces. Using the two lowest-dimensional representations in the auxiliary space, one can eliminate 
the nested Baxter polynomial The resulting equation for the momentum carrying Baxter 

polynomial Q^(u) remains of the second order in finite differences, however, it depends quadrat- 
ically on transfer matrices. The functional algebraic equation in terms of the symmetric transfer 
matrices then reads 

(tW(u + ) - tW(u)tM(u + *))eW" + )Q( 4 )( M + {) (5.1) 
+ (x-) L (tW(u + i) - (x + ) L A^- {u+)+l 2^ {u+) ^j e^ u+)+A -^Q^(u - i) 
- (V 2 V) - (x + )H^(u)e^ u+)+1 2 A ^ u+ ^ ei A -^QW(u) = 0. 
While using antisymmetric transfer matrices, one finds instead 

t [2] {u + )e\^ u+) Q^\u + i) (5.2) 
- (x-) L (t {1] (u + (a; +)^ A -(« + ) + ! A ^ + ^ e ^-(« + H A -(«-)QW( M _ i) 

t [2] {u + )-t {1] {u)t [x] {u + i) + (x + ) L t [1] ( M )e^ A - (M+)+ i A+(u+) ) A A ^Q^\u) = 0. 

As we see, the Baxter equations for the momentum carrying polynomials completely decouple 
from the one built from the nested Bethe roots u k . For the nested Baxter polynomial Q^(u), 
one finds in analogous manner the TQ-relations in terms of t[ a ] 

t [2] (u + )Q^(u + ) = ( x +)i e i A +(« + ) + K(" + ) ^( x+) L e ±A + («+) + iA-(«+) _ f [1](u [+2]^ Q(3)( U ") , 

(5.3) 

and t* s } transfer matrices 

t«(u)t {1} (« [+2] ) - t {2} (« + ))g (3) (n + ) (5.4) 

= (x +)^A + (« + )+lA-(« + ) ( (x+) L e |A + (n +)+ iA_( U+ ) _ t {l }(M[+2]) \ §(8)^-). 



It is straightforward to see that both sets of equations, i.e., Eqs. (I5.ip . (I5.3P and (I5.2p . (15.41) . are 
related to each other by means of the fusion relation for the transfer matrices t\i](u + )t\i](u~) = 
t[2](u) + t^(u). Below in Section El we perform an explicit diagonalization of these Baxter 
equations. 
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5.2 sl(2\2) sector 

Incorporation of both fermionic nodes of the Kac-Dynkin diagram adjacent to its central node, 
with all other being trivial, corresponds to the largest particle-number preserving sector of the 
dilatation operator in the maximally supersymmetric gauge theory. Its fermionic and bosonic 
excitation content consists of two gauginos A 1 , A4, and the complex scalar field X = 42 and the 
covariant light-cone derivative T> + , respectively. 

To write down the Baxter equations it suffices to use three transfer matrices t^(u), t^{u) 
and t^{u) or t[ 2 ](u) and t[ 3 ](u). However, contrary to the rank two sector discussed in 

the previous section, the nesting in the s/(2|2) case is intrinsic and cannot be eliminated. As a 
consequence, the eigenvalues for the momentum carrying Baxter polynomial arise from the 
TQ-relation 

*[i](«) + (^e^*^ 

Q ( - 3) (u+) Q (5) (u-)J 

- (x+) L e^^^Q^\u + ^ (5.5) 

with the dressing factors explicitly depending on the nested Baxter polynomials. In turn, the 
nested Baxter functions obey autonomous Baxter equations, which read in terms of either sym- 
metric 

[t^(u)t^(u^) - tW(u + )tW(u + )] el A + {u+) Q^(u + i) (5.6) 

- (xW) L tW(uW) [*«(«) - tW(u + )tW(u-)} e lA-(^) + |A + ( u [ +2 ])^(3) (w) = 5 
[tW(u)tW(u^)-tW(u-)tW(u-)} e\^ + )Q^( u -t) (5.7) 
+ xV~ 2l ) L * {1} (^ ( V~ 3] )^ A+( ^ 

- (x^)H^(u^) [tW(u) - tW(u + )tW(u-)] el A -^ + l A +^Q^(u) = , 

or antisymmetric transfer matrices 

[t[2](u)t [2] (u^) -t {1] (u + )t [3] (u + )} Q®(u + i) (5.8) 
+ x L (x^)% ] (u + )t [l] (u^^ 

+ (x[ +2 ]) L t [1] ( M [ +3 ])t [2] Hel A -(" [+2, ) + l A +( u[+2] )Q( 3 )H = 0, 

[t[2](u)t [2] (u^) - t {1] (u-)t [3] (u-)} Q ib \u - i) (5.9) 
+ x L {x^)% ] {u-)t [1] {u^)^ 

+ (a;[- 2 ]) L t [1] ( M [- 3 ^ [2 ]Hel A -^ 2] ) + l A+ ^ 2, )g( 5 )( M )=0. 
As we observe the resulting Baxter equations remain of the second order in finite differences. 



10 



5.3 "Gluonic" sector 



Let us also address the aligned-helicity gluonic sector which is closed under renormalization at 
one loop order [12] (see Refs. [59j EO] for recent discussion of related gauge sectors). The space 
of states of the dilatation operator in this case is spanned by single-trace Wilson operators built 
from the twist-one components of the gauge field strength tensor 

F +T = \ (F +1 - iF +2 ) = iF, 2 = -iF u , (5.10) 
F +± ee \ (F +1 + iF +2 ) = -iF\ = iF i{ , 

with either all helicity up F +T or all down F +± and and arbitrary number N of covariant 
derivatives V + , 

tr[£f (F +T ) L ], tr[V*(F +± ) L ]. (5.11) 
In the last equalities in Eq. f)5.10p . we used the spinor decomposition F^ u = -^o^ ' cPFff- + 

^^ ua pF^ 'a which clearly identifies helicity content of the involved Lorentz projections. 

The embedding of this sector into the su(2,2\4) spin chain requires identification of the 
excitation numbers n p associated with it. The use of general excitation numbers from Appendix 
\K\ yields however states which are superconformal descendants of (15.111) . Instead, one can use a 
sequence of particle-hole transformations starting from the distinguished sw(2,2|4) Kac-Dynkin 
diagram in FFFFBBBB grading (with the node numbering from left to right) where the above 
operator corresponds to the excited sixth node with n 6 = N and nontrivial Kac-Dynkin labels 
w 5 = 2 and w 6 = —3 on its fifth and and sixth nodes, respectively. Dualizing the one-loop 
distinguished Kac-Dynkin diagram to the one associated with the long-range spin chain used 
in the present paper, see Fig. [lj one finds the following nonvanishing excitation numbers, = 
2L + N — 2, = 2L — 4 and n-j = L — 2. Of course, beyond leading order of perturbation 
theory, the gluonic operators get mixed with other Wilson operators of different field content and 
possessing identical quantum numbers. 

The construction of the Baxter equations follows the steps outlines in the previous two sec- 
tions: one writes down the lowest-dimensional transfer matrices and solves them, if possible, 
with respect to the nested Baxter polynomials. Eliminating the nested Baxter functions from 
the antisymmetric transfer matrices with lowest three dimensions in the auxiliary space, we find 
the functional relations for the momentum carrying Baxter polynomials Q^(u) which takes the 
form a second-order finite-difference equation 

d a + {u)Q (4 \u + i) + dt(u)Q {i \u -i) + d*{u)Q (4 \u) = , (5.12) 
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where the dressing factors contain residual dependence on the auxiliary polynomial Q^(u), 



d a + (u) = (x + ) L e2 



L jA + (u-)+A + (u+) 



-t [2] (u^) [ (;r [-3] ) L e |A_( U I-3]) + |A + ( U [-3]) + I CT ( B )( u -) 



-(^) L e^ A - ( ^ )+ ^ A+( ^ )+ ^ <5){M[ ~ 3,) 

2hl iA_(«-)4.(5)(«[-3]) 



(5.13) 



dS(«) = e^^") f (x^) £ t w («[- 4 l) [t M («") -f^u^H*!)] e ^-(^) + iA + (^) + i^)(n-) 



- [*[3](^ [ - 2] ) -*[l](^)*[2](^ [ - 3] )] 



^[-3]^ e |A_( u [-3])+|A + ( u [-3])+| CT (5)( u -) 



_( x -)i e i A - (u " )+ i A+(u " )+ i CT(5)(u[ " 31) 



Therefore, this have to be supplemented by TQ-relations for the nested Baxter functions which 
read 

(x + ) L e A+ ^ u+) d a _(u)Q^(u-) - (x~) L e A - {u ~ ) d a + (u)Q {5) (u + ) = , (5.14) 
(x-) L e A - (u ' )+1 2 a(5) ^d a + (u)Q^(u + + (:r+) L e A +( u+ ) + ^ (5) ^cf (u)Q< 6 >(« - i) 

+ e ±A_(«-)-±A + («-)+± CT (5>( U -) h x +)L e A + (u+) d a { ^ +t[1] ^dl( U )\ Q<®(u) = 0, 

Analogously, one finds the Baxter equations for the momentum carrying Baxter polynomials 
Q^(u) in terms of symmetric transfer matrices 

d s + (u)Q (4) (u + i) + d s _(u)Q (4 \u + d s (u)Q (4 \u) = , (5.15) 

while the ones for the nested Baxter functions are 

(x + ) L e A+ ^d s _(u)Q (5 \u-) - (x-) L e A -^d s + (u)Q (5 \u + ) = , (5.16) 

. -)L e AAu-)+^(u+) d s^ u)Q (6)^ u + t) + ( x +)^A + (« + )+|^(«-)^ (M) g(6) (M _^ 

+ e |A_(«-)-iA + («-)+iaW(«-) r (a . + ) L e A + ( U +) rf s (M) + t W (u)d 8 + (uj\ Q (6) («) =0, 



[X 
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with 

d s Ju) 



d s _{u) 



d s (u) 



yX J e 2 



(x-) i ^ 1 >(n[- 2 I)^ 1 >(«[- 4 ])e2 A -( w -)+I A +( w -)+2- (6) (« [ - 3] ) 



(x [ ~ 3l ) L e^ A - ( " [_31)+ ^ A+( " [ " 3,)+ ^ (5)(u " ) 



( x -^ e |A_(n-)+iA + («-)+i CT (5)( u [-3]) 



[t{ 2 >( U -)tW( U [-3]) _t{l}( u [-2]) t {3}( w [-2])] e , 
1 

-e2 



2]jj e |A_(u-)+|<7C 6 )(«I- 3 ]) ^ 

§ A +^ | (x-) i ^ 1 >( w [- 4 i)^ 2 >( w -)ei A -( u ") + l A +^) + i CT(6) ^ 3] ) 

(x h3l ) L ei A - (n[ " 3,)+ i A+( " [ " 31)+ ^' j(5)(M " ) 

-(x") L ei A - (u ~ )+ l A+(n_)+ i CT(5)( " [ " 31) 



(5.17) 



-t< 3 >( u [- 2 ]; 



As anticipated, for vanishing coupling constant, g — 0, the nesting disappears indicating that 
the aligned-helicity gluon operators form a closed subsector of the dilatation operator 



6 Solution of fusion hierarchies 

Finally, let us address the solution of the fusion hierarchy for the rank-two s/(2|l) sector, i.e., 
finding explicit eigenvalues for all transfer matrices for a given state of the spin chain. 

The infinite tower of transfer matrices t\ a ] (u) and t^(u) can be expressed in terms of just two 
low-dimensional antisymmetric transfer matrices £rii(w) and t[2](u). Namely, using the analogue 
of the duality relation (13.71) for the sZ(2|l) sector and solving it with respect to t[s\(u) as in Eq. 
(14.71) . one finds 

, / \ t [2] {u-)t [2] {u + ) 

tmw = 1 — — — — - (6-i) 

t m {u) - (x-) L e2 A +(^) + 2 A -(^) 



Q(3)( 



II. 



QW(ul 



-i 



All other transfer matrices with antisymmetric auxiliary space can be found iteratively 

= ■^ ) ^t^^) l W ^ , (6 . 2) 

- - Q(3)( M [-a+l]) 

for a > 3. While all transfer matrices with symmetric representations in the auxiliary space 
follow from the Bazhanov-Reshetikhin formula (13.41) . 

Therefore, the problem is reduced to the solution of functional relations (15.21) and (15.31) for 
Baxter polynomials Q^\u) and Q^ 3 '(u) and transfer matrices tu](u) and t\2](u). At leading order 
of perturbations theory, i.e., for g = 0, Eq. (15.21) can be solved immediately along the same lines 
as the si (2) Baxter equation analyzed in Appendix O However, beyond one loop its successful 
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solution requires proper identification of the renormalized spectral parameter which serves as 
an expansion variable for polynomial fused transfer matrices. The reason for this is that the 
perturbative series for the transfer matrix t\ a ] acquires nonpolynomial terms in the bare spectral 
parameter u. This happens before the wrapping order 0(g aL ) sets in, when the asymptotic 
Baxter equation is expected to work and this is ensured by the right choice of the expansion 
variable. A guidance into the acceptable choice can be gained by comparing Baxter equations 
in terms of fused transfer matrices with the ones written by means of conjugate fundamental 
transfer matrices introduced in Ref. [35]. Namely, in terms of the transfer matrix £m and its 
conjugate tji], one finds instead of Eq. (15. 3p 

^Ja+(«-)+!a_(«-) _+„f„^\ n(3)/,. + \ _ ( / p |a+(«+)+±A_( U +) _Tt.,\ \ n(3)(„-^ 



(x-) L e2* +(u > + 2*-"> >-t {1] (u)jQW(u+)= I (x + ) h e2*^ >+ 2*-^> -i {1] (u)) Q^(u 

(6.3) 

Eliminating from Eqs. (I5.3P and (16.31) . one deduces 

U u) = ^eW-Hi^M ^ei^HK^) _ ^ Z ^) 

V / xL e 2 A +M+2 A -M -t m (u-) 

(6-4) 

which suggests that the transfer matrix t[2](x) can be factorized into a product of two functions, 
one depending on x and another on x' - ^. For states with real charges 0^ entering t\x](u), 
3m [0^] = 0, it further simplifies to 

t [2] (u) = a^+M+KM ^ x [- 2 ]^ e |A + («[-]) + iA_(.M]) _ 

Taking into account these considerations, the transfer matrices admit the following expansions 
in the renormalized spectral parameters 

«[!](«) = £ alto (e-s) 

fc>0 

t M («) = J2* L ~ j *j ] (9) J2^ ] (9) {x [ - 2] ) L ~\ (6.6) 

j>0 k>2 

where the tower of JH!- —charges encodes the expansion of the exponential dressing factors in Eq. 
(16.41) . The first few charges in both transfer matrices (I6.5P and (I6.6P can be found in a closed 
form in terms of the Baxter polynomial Q^\u) by studying the asymptotics of tm(u) and tm(u) 
via the fusion relations (12.111) . This requires the large— u behavior of the dressing factors A±(u). 
Using Eqs. (12.61) and (12. 9p . one finds that the first few terms in the 1/u— expansion yields 

u u z 

where the residues are (a = 1,2) 



A H = y ^y/T-filw^{fft,g) (lng( 4 )(±f - - #W(fft,g) fin 



Q (4) (-f-^) 



(6.{ 
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and the explicit form of the functions entering the integrand is 



(t, g) = -g 2 , (t, g ) = , (6 .9) 

g 2 



w^(t,g) = y 2 t, #W(t,g) = 8iZ 2 , 1 (g) 

with |M EH 



From these considerations, one immediately concludes that Q^(g) = (gr) = 1 for both a = 1, 2 
and the all-order conserved charges 0^ and 0^ °f the fundamental transfer matrix tm(w) read 

0[ 11 (^) = |(AV 1) (^) + A (1) (^)) , (6.11) 

£% ] (g) = -Q®(g) + § (0' 1] (S)) 2 + \ (a^) + A®(gj) , (6.12) 

[21 

with the latter charge expressed in terms of 0[ introduced below. On the other hand, for the 
renormahzed fused transfer matrix Tm{u) = x~ L t[ 2 ]{u), one deduces by the same token that the 
leading asymptotics as u — > 00 

T m {u) = C 2 (g)u L - 2 + 0(u L ' 3 ) , (6.13) 
is driven by the value of the renormahzed quadratic Casimir of the s/(2|l) algebra 

C 2 (g) = (n 4 + §7(0)) (n 4 - n 3 + L - 1 + § 7 (#)) , (6.14) 
where the anomalous dimension (I2.14p is expressed in terms of residues of the dressing factors 

7 (^) = ^(A ( h 1) (^)-A (1) (^)) . (6.15) 

This charge, as we saw above, defines up to an additive contribution, the subleading term in 
the fundamental transfer matrix. With this knowledge at hand, we can understand why the 
expression in braces in Eq. (16 Ah does not have the induced charge 0^: it cancels between the 
two terms in the large— u limit, however, it re- appears in the multiplicative factor in front of the 
polynomial in a^ -2 ' variable, see Eq. (16.61) . So that 

Qg ] (g) = C 2 (g), X l ?(g) = £l®(g). (6.16) 

This information suffices to solve the multi-loop Baxter equation to arbitrar)Q order of per- 
turbation theory. In Table [1] we give, as an example, the lowest four orders in the expansion of 
the anomalous dimensions and conserved charges 

iM = 5> 2£+ v,, 0L 1] (s) = E^il (6-i7) 

e=o £=0 

for the state a = [L, 713,714] = [5,3,7] of the sl(2\l) long-range spin chain. Conserved charges 
associated with other fused transfer matrices can be found by means of Eqs. (16 .4p . (16 .ip and 



4 We implicitly imply that it holds only below the wrapping order due to the asymptotic character of Baxter 
equations. 
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k 





1 


2 


3 




7 


679 
96 


59297 
4608 


6355979 , 175 (■ (o\ 
221184 ' 96 > v / 


19 

3 


5587 
864 


735545 
62208 


473456465 , 493 A /q\ 
17915904 t 288 ^ v ^ 







12 


83^ 
1152 


27703^ , V21s/o\ 
221184 48 °) 





V7 
12 


1025-/7 
12096 


381 1897^7 V7/-(o\ 
24385536 144^°/ 




-56 
-56 


105 
2 

95 
2 


1295 
32 

229 


36267 7fiQ\ 
512 

465511 19 / ^'("?^ 
6912 24>V°J 




-8y/21 


79v / 2T 
6 


12887v / 21 
2304 


556337^21 351/21/^0^ 
55296 24 V > 


-16-s/7 


845^ 
42 


268805^ 
21168 


1862658941^7 25 v / 7a/o\ 
85349376 72 ^\°> 




70 

38 


329 
2 

547 
6 


2933 
192 

11401 
1728 


528031 1 63^^ 
9216 + 4 

3761615 | 433a/qN 
124416 36 ^ v°/ 


oS\ 9 ) 


8^21 



6 



19865^ 
2304 




61963^ | 7v / 2Ta ( 'cA 
55296 1 2 





Table 1: Eigenvalues of the a. = [5, 3, 7] states up to four-loop order in Af = 4 super- Yang-Mills 
theory. 
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7 Conclusions 



In this work, we have suggested the hierarchy of fused transfer matrices for the putative long- 
range magnet of the maximally supersymmetric gauge theory. These allow one to construct a 
finite set of Baxter equations for eigenvalues of Baxter polynomials for momentum carrying and 
auxiliary Bethe roots. We have demonstrated the formalism using low-rank subsectors of the 
dilatation operator. 

The T-system of transfer matrices formulated in Section [3] is a starting point for the derivation 
of the Y-systems [371 1381 [39] which, in turn, yield equations of the Thermodynamic Bethe Ansatz 
[35l|36]. The latter, implemented for the mirror image of the long-range spin chain, incorporate 
finite-size effects [631 [M] with the lowest order correction reproducing the Luscher formula (see, 
e.g., Ref. [62]), which was recently reincarnated within the gauge/string duality in Refs. [63] to 
account for wrapping effects evading the framework of the asymptotic Bethe Ansatz. Therefore, 
one anticipates that the framework outlined here may be used to consistently incorporate this 
kind of finite-size effects. 

This work was supported by the U.S. National Science Foundation under grant no. PHY-0456520. 
We are indebted to Romuald Janik, Gregory Korchemsky, Vladimir Korepin, Rafael Nepomechie, 
Adam Rej, Matthias Staudacher, Zengo Tsuboi and Stefan Zieme for useful discussions and 
correspondence at different stages of this work. 

A Superconformal algebra 

Let us present the Serre-Chevalley basis associated with the Kac-Dynkin diagram in Fig. [TJ 
As a corollary of our considerations we will identify the excitation numbers defining Baxter 
polynomials with quantum numbers of gauge invariant Wilson operators that they correspond 
to. We start with the su(2, 2|4) superconformal algebra in the notation of Ref. [16] and use spinor 
projections 



The chirality charge R is merely identified R = 01. The inverse transformations for the bosonic 
generators of the so(4, 2) conformal algebra are 



= -\ tr [a^ L] + \ tr [a^ L] , = tr [a* P] , i& = tr [a* K] . (A.2) 



QaA — ^O-aA , Qf = ±0f, 



^ /3 ' 



9£ = i6i, 




The generators (lA.lj) can be organized in a supermatrix 



/ ^°/3 2 



(D + iR) 



\ 



Q 



■0 



QaB 



^A 



(A.3) 



V 




hj + \8J (D - |R) J 
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with indices running over eight values A 
relations 



(a, A, at). They obey the graded commutation 



[E A B ,E e v} 



E A B E C V 
S C B E\ - 



>_ 1 yA+B)(c+v) E c vM A B 
_ l) {A+BW + f>) 5 A vK c B _ 



(A.4) 



The root system of the su(2, 2|4) algebra is expressed in terms of weights which form a basis 
of the dual Cartan subalgebra. We enumerate them as follows 



v a = (e 3 , e 4 |<5i, $2, 8 3 , <5 4 |ei, e 



and endow them with inner products 

5, 



i £i £ i 



Hj 1 



(8i\8j 



Hj 1 



2) , 



0. 



The simple root system corresponding to the Kac-Dynkin diagram in Fig. [T] reads 



OLl 



= £2 
OC5 ~- 



-8l, OL 2 

£3 - 8 3 , 



81-82, OL 3 

olq = 8 3 - 84 



and the resulting Cartan matrix A pq 

( 



A 



(OL p \OL q 
1 

-1 



8 2 -£1, 0:4 = 6! 
CK7 = 84 — £ 4 , 



\ 



£3, 



V 



1 

-2 
1 



defines the Serre-Chevalley basis of the superconformal algebra 



hi 


= E 8 8 + 


E 3 3 = 


-Li 1 + iB — ly , 


4 


= E 8 3 = 


Q21, 


e r 


= E 3 8 = 


§ 21 , 


h 2 


= -E 3 3 


+ E 4 4 


= Tx 1 -T 2 2 , 


4 


= E 3 4 = 


-T 2 \ 


e 2 


= E 4 3 = 


-T1 2 


h 3 


= -E 4 4 


-E 7 7 


= -L 1 1 -i© + T 2 2 , 


4 


= E 4 7 = 


§ 12 , 


e 3 


= E 7 4 = 


Q 12 , 


hi 


= E 7 7 - 


E 1 ! = 


Lx 1 + L\ +D, 


4 


= E 7 x = 




e 4 


= E 4 7 = 


— K 41 


h 5 


= E\ + 


E 5 5 = 




4 


= E 4 5 = 




e 5 


= E 5 i = 


Of, 


h 6 


= -E 5 5 


+ E 6 6 


= T 3 3 -T 4 4 , 


4 


= E 5 6 = 


-T 4 3 , 


e 6 


= E 6 5 = 


-T3 4 


h 7 


= E 2 2 - 


E 6 6 = 


- ±B + T 4 4 , 


4 


= E 6 2 = 




e 7 


= E 2 6 = 


-si, 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



via the conventional commutation relations [h p , = AzA pq e^. Notice that other simple root 
systems are related to this one by reflections with respect to the fermionic simple roots cx p 
(p = 1, 3, 5, 7) with the vanishing bilinear form (ex lex) = 0. They span the Weyl supergroup. For 
Bethe Ansatz equations these correspond to particle-hole transformations. 

Making use of graded commutation relations (1A.4I) . one can cast the generators (lA.lj) into the 
oscillator representation [65]. The latter is built from the bosonic (a a , a*, b@, &t) and fermionic 
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(ca, c* a ) creation/annihilation operators and reads 

K$ = °&1> K^ = a?b°, B = ±a\a?> + ±b\tf + 1 , (A.10) 

L/ = a\pP - \^a\o? , L* $ = 6*6* - ^b\b"< , 

QaA = c A ai , Q A = c tA 6t , S aA = J A a a , ^ = c A b« , 

T A B = c A J B - \5 B A c D J D , R = 6* - a* a* , 

with the vanishing central charge condition c JfD CD+al / a' r —bl r b A/ = 2. This representation allows one 
to easily identify the excitation numbers n p with eigenvalues of Cartan generators. Introducing 
a generic state of the spin chain, 

\O m ,..,n 7 ) = ( fl ( e 4 +fc ) n4+fc (e 4 + _ fc ) n4 - fc j (e 4 + )" 4 |fi) (A.ll) 

)«5-n6( c tl)n 2 -n^ c t2^n3-n 2 ^ a tjn4-n3( a t)ni( & t)n4-n s ^t^n 7 |^ ^ 

where we remind that the vacuum state is |Q) = (c i " 3 c t4 ) jL |0) = tr[Z L (0)], one observes that 
the excitation numbers n p can be "counted" with the Cartan generators by relating them to the 
number-of-particle operators 




a \a} = i© + Li 1 - \R - \ , a|a 2 = §B - Li 1 - ±R - \ , 

b\b l = lB + L\ + \R-l, b\tf = iD-L i i 



2 1 4 2 

fjj.-i - !-Fn fl. iB> 1 7>"i"7,2 — 1 TTfl f~ 1 . ± ||_ i 

c ti Cl = -T/ + ±R + \ , C t 2 c 2 = -T 2 2 + iR + I , 

CsC t3 = Ts 3 _ 1 R + I ? C4C f4 = T 4 4 — |R + I , 

r4 



(A.12) 



where J2 i=1 TV = 0. This yields the following identities 



ni = 


¥ 


— s — \r—\L, 




n 2 = 




-s-h, 




n 3 = 


¥ 


-s + \r-t 1 -t 2 + 


¥, 


n 4 = 


d- 


■h-t 2 , 




n 5 = 


¥ 


-s-\r + h + t A + 


1 T 


n 6 = 


¥ 


- s + 1 4 , 




n 7 = 


¥ 


-s + \r-\L, 





(A.13) 



where d, s, s, r and U are eigenvalues of D, Li 1 , L 1 ^ R and TV, respectively. These numbers 
obey obvious inequalities, ni < n 2 < n 3 < n 4 > ns > n 6 > 717, L > n 6 — n 7 and L > n$ — 11$. 



B Young tableaux for FBBFFBBF grading 

Considerations presented in the main body of the paper can be applied in a straightforward 
manner to the Kac-Dynkin diagram in the FBBFFBBF grading, i.e., with the compact central 
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su(2) node. The elementary Young tableaux read in this case 



X 



X 



X 



X 



X 



X 



X 



X 



Q {1) (u 


h ) e IK 4) («-)- 


Q (1) (u- 


-) 


Q {1) (u 




{u- 


-) Q ( ' {u) 


Q (2) (u 


+ i) (u~) 


Q { ) ( 




Q (3) (u 


-) g( 4 ) (u + i) 


g( 3 ) (u- 


v) Qw (u) 




h ) (w - i) 


(u- 


-) <2 l ' {u) 


Q [b > (u 


- i) Q (5) (u + ) 


QW ( 




g (7) (u- 


") (u + i) 


(u- 




Q {1) (u 


-) e |p (M+) - 


(ir 





(B.l) 



(B.2) 



The grading of each node is inverted compared to the Kac-Dynkin diagram in Fig. [TJ i.e., 1 = 
4 = 5 = 8 = 1 and 2 = 3 = 6 = 7 = 0. The fused transfer matrices are built along the procedures 
outlined in Sections [2] and [3J 



C Perturbative solution of sl(2) Baxter equation 

For the reader's convenience, we present here a sample Mathematica routine to compute anoma- 
lous dimensions in the si (2) sector to two loop order from the Baxter equation 

(x + ) L e A + ( " +) Q (4) (« + •) + (arjVMtOgW^ _ ,) = t (u)Q w (u) . (C.l) 
The perturbative expansion for Baxter polynomials and transfer matrix 

Q^\u) = Q$\u) + g 2 Q?(u) + 0(g 4 ) , t(u) = t (u) + g%(u) + 0(g 4 ) , (C.2) 
and corresponding expansion of the dressing factors is implemented in the following module: 

Spectrum [spin_, length.] := Module[{J, L, BaxEqO, BaxEql, k, n, Q40, Q41, QQO, QQ1}, 
J = spin; L = length; 

Q40[u_] = u~J + Sum[a[k] u~k, {k, 0, J - 1}]; 
Q41[u_] = Sum[b[k] u"k, {k, 0, J - 1}]; 

t0[u_] = 2 u~L - ((J + L/2)(J + L/2 - 1) + L/4 ) u~(L - 2) + Sum[q0[k] u" (L - k) , {k, 3, L>] ; 
tl[u_] = Sum[ql[k] u"(L - k) , {k, 1, L>] ; 

BaxEqO := Q40 [u] t0[u] - (u + 1/2) ~L Q40 [u + I] - (u - 1/2) ~L Q40 [u - I]; 
BaxEql := (Q41 [u] t0[u] + Q40 [u] tl[u]) 
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- (- (L/4)(u + 1/2) "(L - 2) Q40[u + I] + (u + 1/2) "L Q41 [u + I] 

- (u + 1/2) "(L - 1) Q40[u + I] Q40'[I/2]/(2 Q40[I/2])) 

- (- (L/4)(u - 1/2) "(L - 2) Q40[u - I] + (u - 1/2) "L Q41 [u - I] 

- (u - 1/2) "(L - 1) Q40[u - I] Q40' [-1/2]/ (2 Q40 [-1/2] ) ) ; 
QQO = Q40[u] ; QQ1 = Q41 [u] ; 

BaxEqO = Collect [BaxEqO , u] ; BaxEql = Collect [BaxEql , u] ; 

Do[SolO = Solve [Coefficient [BaxEqO, u, L + J - 3 - (k - 1)] == 0, a[J - k] ] ; 

BaxEqO = Simplify [BaxEqO/ . SolO [ [1] ]] ; QQO = Simplify [QQO/ . SolO [ [1] ]] ; 
BaxEql = Simplify [BaxEql/ . SolO [ [1] ]] , {k, 1, J}]; 
SolutionQQex = NSolve [Table [Coefficient [BaxEqO , u, k] == 0, {k, 0, L - 3}], 

Table[q0[k + 3], {k, 0, L - 3>] , 100]; 
Do[SolQll = Solve [Coefficient [BaxEql, u, L + J - k] == 0, ql [k] ] ; 
BaxEql = Simplify [BaxEql/ . SolQll [ [1] ]] ; 
QQ1 = Simplify[QQl/.SolQll[[l]]] , {k, 1, L}] ; 
Do[Soll = Solve [Coefficient [BaxEql , u, J - k] == 0, b[J - k] ] ; 
BaxEql = Simplify [BaxEql /. Soll[[l]]]; 
QQ1 = Simplify[QQl/.Soll[[l]]] , {k, 1, J}]; 
EnO = Simplify [Factor [(I* (D [QQO, u]/(2 QQ0))/.u -> 1/2) 

- (I*(D[QQ0, u]/(2 QQ0))/.u -> -1/2) ]]; 
Enl = Simplify [Factor [(I*(D[QQ0, u]~3/(8 QQ0~3) 

- 3 D[QQ0, u] D[QQ0, {u, 2>]/(16 QQ0~2) 

+ D[QQ0, {u, 3}]/(16 QQO) + D [QQ1 , u]/(2 QQO) 

- D[QQ0, u] QQl/(2 QQ0~2))/.u -> 1/2) 

- (I*(D[QQ0, u]~3/(8 QQ0~3) 

- 3 D[QQ0, u] D[QQ0, {u, 2}]/(16 QQ0~2) 

+ D[QQ0, {u, 3}]/(16 QQO) + D [QQ1 , u]/(2 QQO) 

- D[QQ0, u] QQ1/C2 QQ0~2))/.u -> -1/2) ]]; 
LogThetaO = Log [Factor [(QQO /. u -> I/2)/(QQ0/.u -> -1/2)]]; 

Dat = {{GammaO -> Rationalize [ EnO/ . SolutionQQex [ [1] ]] , 
Gammal -> Rationalize [Enl/ . SolutionQQex [ [1] ]] , 
logTheta -> Chop [N [LogThetaO/ . SolutionQQex [ [1] ]]]» ; 
Do [Dat = Join [Dat, {-[GammaO -> Rationalize [EnO/ . SolutionQQex [ [n] ]] , 
Gammal -> Rationalize [Enl/ . SolutionQQex [ [n] ]] , 
logTheta -> Chop [N [LogThetaO/ . SolutionQQex [ [n] ] ] ] »] , 
{n, 2, Length [SolutionQQex] }] ; 

Dat] 

The output is presented in the form {GammaO— ■> 7 , Gammal— > 71, logTheta— > ln#}, with 70 and 
71 being the numerical values of anomalous dimensions at one and two loops, respectively, and 
log^ is the value of the corresponding quasimomentum 9 of the state. Though we used the 
explicit form of the conserved charge q0 [2] in terms of eigenvalues of the bare quadratic Casimir 
C 2 (g = 0) = -(N + \L)(N + \L - 1) - \L of the sl(2) algebra where n 4 = N is the number 
of excitations and L is the length of chain, cf. Eq. f)6.14p . it is not a necessary input as its value 
can be determined from the Baxter equation itself as a coefficient in front of the highest power 
of the spectral parameter u L+N ~ 2 . Generalization to higher loop orders is straightforward and 
analyses were presented elsewhere [44] . Baxter equations for higher rank algebras can be treated 
analogously with results for s/(2|l) long-range spin chain presented in Section El 
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